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We present observations of two poor clusters, AWM 4 and MKW 4, observed by 
XMM-Newton. Both systems are relaxed, with little substructure evident in their X- 
ray halos or galaxy populations. However, their temperature structures are markedly 
different, with AWM 4 isothermal to the resolution of the EPIC instruments while 
MKW 4 shows a strong decline in temperature towards the core. Metal abundance 
also increases more strongly in the core of MKW 4 than AWM 4. Three dimensional 
models show further differences, suggesting that gas in the core of AWM 4 has been 
heated and has expanded outwards. The dominant elliptical galaxy of AWM 4 hosts 
an AGN with large-scale radio jets, while the central cD of MKW 4 shows no current 
AGN activity. We therefore conclude that the difference in activity cycle of the AGN 
in the two galaxies is responsible for the difference in IGM properties between the 
two clusters. 



1. Introduction 

MKW 4 and AWM 4 are fair ly relaxed poo r clus- 
ters, originally identifie d in the iMorgan et alJ l(l975l ) 
and lAlhert et all l|1 977ft surveys. Each is dominated 
by a giant elliptical or cD galaxy, surrounded by —50 
(MKW 4) or -30 (AWM 4) other galaxies. Kinematic 
studies of the galaxy populations of the two clusters 
show them to be fairly evenly distributed about their 
mean redshifts, with no sign of significant substructure 
llKoranvi fc Gellerll2002ft . Both clusters are morphologi- 
cally segregated, with absorption line systems found pre- 
dominantly toward the center of each system. 

The dominant galaxies of MKW 4 and AWM 4 are 
NGC 4073 and NGC 6051 respectively. Both are highly 
luminous and ~1.5 magnitudes brighter than the second 
ranked galaxy in their clusters. NGC 4073 shows evi- 
dence of_a_p_as^nBrg_er. It has a kinematically decoupled 
core l)Fisher et al.lll995l) and t he spectroscopic age of th e 
stellar population is —7.5 Gyr l)Terlevich &: ForbeJ2 002'). 
It is the larger of the two galaxies, with Log Ls = 11.01 
Lq. NGC 6051 hosts an active nucleus, revealed by a 
large scale radio feature extending from the galaxy core. 
The central point source is not detected at other wave- 
lengths, which suggests that the axis of the AGN jets 
is in the plane of the sky and that the central engine is 
highly absorbed. The optical luminosity of NGC 6051 is 
£5=10.76 2,0. 

Both clusters have been the subject of previous X- 
ray observations with Einstein, ASCA and ROSAT, 
all of which have shown them to be relaxed sy stems 
ijdell'Antonio et alJ Il995t Uones fc Formanl 1 19.991) rea- 
sonably well modeled by at most two beta models 
ijHelsdon fc Ponmanll200a iFinoguenov et al.ll200lfl . The 
X-ray halo of MKW 4 has a mean temperature of —1.7 
keV, and radial temperature profiles show that the tem- 
perature increases from —1.3 keV in the core to a peak 
above 2.5 keV at —115 kpc before dropping back to a 



fairly constant temperature of —1.5 keV at higher radii. 
The ASCA data have been used to produce radial abun- 
dance profiles which show Fe, Si and S abundances which 
fall from central values of 0.4-0.7 to 0- 0.2 Z Q at a ra- 
dius of 500 kpc l(Finoguenov et all 120001. The halo of 
AWM 4 is somewhat simpler, having a relatively isother- 
mal temperature profile out to a radius of at least —200 
kpc. The mean temperature of the gas is —2.5 keV. 

Given their relatively relaxed and undisturbed state, 
and the large mass concentrations associated with their 
central galaxies, these clusters are potentially excellent 
sites in which to observe cooling at the low end of the 
cluster mass range. We have used two XMM-Newton ob- 
servations to study the structure of their X-ray halos, and 
in particular to determine whether they currently host 
cooling flows. Throughout these proceedings we assume 
Hq=75 kms^Mpc^ 1 and normalise optical luminosities 
to the B-band luminosity of the sun Ls0=5.2xlO 32 erg 
s _1 . MKW 4 lies at a redshift of z=0.020, AWM 4 at 
2=0.032, and we therefore assume distances to the two 
clusters of 80 Mpc and 127 Mpc respectively. 

2. Results 

MKW 4 was observed with XMM-Newton during or- 
bit 373 (2001 December 21) in two exposures of -16,000 
and —4,500 seconds. AWM 4 was observed during orbit 
573 (2003 January 25-26) for just over 20,000 seconds. 
In both cases the EPIC MOS and PN instruments were 
operated in full frame and extended full frame modes re- 
spectively, with the medium filter. A detailed summary 
of the XM M-Newton mission a nd instrumentation can 
be found in lJansen et al. I fcOOll) . and references therein. 
The raw data from the EPIC instruments for the longer 
of the two exposures were processed with the publicly re- 
leased version of the XMM-Newton Science Analysis Sys- 
tem (SAS v.5.3.3 for MKW 4, SAS v. 5.4.1 for AWM 4), 
using the EPCHAIN and EMCHAIN tasks. After filtering 
for bad pixels and columns, X-ray events corresponding 
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to patterns 0-12 for the two MOS cameras and patterns 
0-4 for the PN camera were accepted. Investigation of 
the total count rate for the field revealed a short back- 
ground flare at the beginning of the MKW 4 observation 
and another in the second half of the AWM 4 observa- 
tion. Times when the total count rate deviated above 
the mean by more than 3<r were therefore excluded. The 
effective exposure times for the MOS and PN cameras 
were 14.1 and 10.5 ksec in the case of MKW 4 and 17.4 
and 12.7 ksec in the case of AWM 4. Images and spec- 
tra were extracted from the cleaned events lists with the 
SAS task evselect. FigureQshows adaptively smoothed 
images of the two clusters. 

2.1. Surface brightness profiles 

In order to model the surface brightness distributions 
of the two clusters, we produced images for use in the 
CIAO SHERPA fitting software. Using data from all three 
EPIC cameras we searched for and removed all point 
sources, with the exception of false identifications associ- 
ated with the cluster cores. We also generated back- 
groun d im ages based on the blank - sky data of iLumbl 
( 2002) and Plead fc PonmarJ ll2003ft . and the telescope 
closed data of iMartv et alJ ll2002i). using the 'double 
subtr action' technique l)Arnaud et al J l2002i iPratt et al.l 
l200ll) . As the images contain many pixe ls with few (or 
zero) counts, we used the Cash statistic l|Cashlll97f^ to 
perform the fits. 

AWM 4 was reasonably well fit by two beta models, a 
more extended elliptical component and a circular cen- 
tral component. The parameters (and la errors) of the 
best fit model are shown in Table ^ Errors for each 
model were estimated with all fitted parameters free for 
perturbation. Position angle is measured anti-clockwisc 
from north. 

We initially attempted to model MKW 4 using either a 
single or two-component beta model, but this was unsuc- 
cessful. We decided instead to model sectors of the halo 
individually. We selected four sectors, covering the ma- 



jor and minor axes of the halo. The NW and SE sectors 
showed deviations from a simple beta model, with the 
NW sector showing emission extending to higher radii, 
and the SE sector showing a 'bulge' of excess emission 
at ^20 kpc. The two remaining sectors were well de- 
scribed by a single beta model and had similar core radii 
and slopes, leading us to fit them together. Taking into 
account the effect of the PSF and excluding all data out- 
side the NE and SW sectors, we found the best fit to 
be a beta model with parameters as shown in Table ^ 
Note that this model is by necessity circular, with the 
axis ratio fixed at 1.0. 

2.2. Temperature and Abundance profiles 

In order to examine the large scale temperature and 
abundance structure of the clusters we produced radi- 
ally binned spectra of the halo. Circular bins were used 
for MKW 4, elliptical bins with ellipticity and position 
angle taken from the surface brightness fits for AWM 4. 
All data within 17" of point sources were removed and 
spectra were fit using an absorbed MEKAL model in a 
0.4-4 keV energy band. Figure shows the results of our 
fits. 

The deprojected temperature profile for MKW 4 shows 
a fairly smooth decline in temperature towards the core, 
consistent with ROSAT results (Helsdon & Ponman 
2000). Metal abundance is strongly affected by changes 
in hydrogen column. The fit quality in the central bins 
is poor, and we experimented with different models to fit 
these bins. The fits could be improved by using a two- 
temperature model, but the higher temperature compo- 
nent was never well constrained by the data. Using cool- 
ing flow models did not produce good fits except in cases 
where the emission from low temperature gas (i.e. gas 
at temperatures lower than ~ 1 keV) was minimal. 

The deprojected temperature profile for AWM 4 shows 
a fairly constant temperature, with a possible slight de- 
cline with increasing radius. To test whether this decline 
is real, we tied the four inner and two outer bins to form 
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Fig. 2. — Deprojected temperature and abundance profiles of the two clusters. Symbols represent 90% error regions. Fitting was carried 
out using the XSPEC PROJCT model, and assuming the emission arises from an absorbed MEKAL model. Hydrogen column was fixed at 
the galactic value. The grey regions in the AWM 4 plots (upper panels) show fits with bins tied together to form inner and outer region. 
Note that the temperature difference is still not significant. 



two large spectral regions. The results from these are 
marked as grey boxes, and show that the decrease is not 
statistically significant. The abundance profile shows a 
more marked decrease with radius. The spectra from all 
bins were well described by single temperature plasma 
models, and no more complex models were required. 
For each cluster, we also investigated the abundances 



of certain individual elements, using projected spectra so 
as to avoid reducing the quality of the data being fitted. 
Using the VMEKAL model in XSPEC we are able to fit 
Fe and Si (and in the case of MKW 4, S) individually, as 
well as the abundance of remaining metals and the tem- 
perature in each bin. Figure y| shows abundance profiles 
for Fe, Si and S in MKW 4, shown with reference to the 
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Fig. 3. — Radial abundance profiles for Fe, Si and S in MKW 4. 
In each plot, diamonds show the 90% error on the metal in question, 
as labeled on the plot. The crossed error bars show the fit to all 
the remaining metals, again with 90% errors. Njy was fixed at the 
galactic value in all fits. 
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Fig. 4. — Radial abundance profiles for Fe and Si in AWM 4. 
Diamonds show the 90% error on Si, Boxes represent Fe and the 
crossed error bars show the fit to all the remaining metals, again 
with 90% errors. Nh was fixed at the galactic value in all fits. 



profile for all other metals combined. There is a clear 
increase in all three elements in the core of the cluster, 
at the location of NGC 4073. Figure 0] shows abundance 
profiles for Fe and Si in AWM 4, and again shows an 
increase in abundance toward the core, at the location of 
NGC 6051. 

Based on the elemental abundances of Fe and Si found 
in the inner bins of the two profiles, we are able to es- 
timate the contribution made by SNIa and SNII to the 
enrichment of the gas in the cores of the two clusters. 
We assume yields of Si and Fe for t ype II supernovae 
of 37^ =0-07 M and y Si =0.133 M Q l)Finoguenov et all 
2000), and from type la supernovae ;/fy =0.744 M Q and 
2/ Sl =0.158 M Q ^Thielemann et alj fl993). For MKW 4 
we estimate that ^50% of the iron within 17.5 kpc of 



the core was injected by SNIa. This is similar to our 
result for AWM 4, where we find that ~52% of the iron 
within 49.2 kpc was injected by SNIa. However, we note 
that these values are strongly affected by our choice of 
definition of solar abundan ce ratios. We use the ratios of 
Anders fc Grevessel |l989jl , but the more recent ratios of 
Grevesse fe SauvaJ l)l998|) give a rather lower value for 
Fe/H (Fe/H=3.2xl0" 5 compared to 4.7xl0~ 5 ). This 
would suggest that our Fe abundances are underesti- 
mated by a factor of ~1.4 compared to those calculated 
using the more recent ratios, and similarly that our esti- 
mates of SNIa contribution are lower than those we would 
find were we to use the Grevesse & Sauval abundances. 

2.3. Three dimensional gas models 

Using software provided by S. Helsdon, we use the sur- 
face brightness fits derived in Section l2~D and the radial 
temperature profiles shown in Figure|21to calculate quan- 
tities such as gas mass, total mass, gas cooling time and 
gas entropy for the inner region of each cluster. The 
software infers the gas density profile based on the fit- 
ted surface brightness profile and a model description of 
the temperature profile, and normalises it by comparison 
with the X-ray luminosity. Given this density profile we 
can use the well known equation for hydrostatic equi- 
librium to calculate the total mass at any given radius. 

2 

Entropy is defined as S = T/n| . Figure [5] shows radial 
profiles for these and other relevant parameters for both 
clusters. 

Although the two clusters have quite different struc- 
tures, their total mass profiles are very similar outside the 
core regions. This implies that the two clusters have sim- 
ilar total masses, with the difference in the core probably 
caused by the different masses of the dominant galaxies. 
Assuming a mass-to- light ratio for stars of 5 M /L®, we 
would expect the two galaxies to have total stellar masses 
of ~5xlO n M Q and -3xlO n M , values which are com- 
parable to the total mass of each system at 10 kpc. 
MKW 4 has a larger total mass and a larger mass of 
gas in its core regions than AWM 4. The gas fraction 
of AWM 4 is very low out to ~40 kpc, where it begins 
to increase rapidly. In the outer regions the situation 
is reversed, and AWM 4 has a considerably larger gas 
fraction than MKW 4. MKW 4 has a relatively short 
cooling time in the core (~5xl0 8 yr), and the inner halo 
contains very low entropy gas. In AWM 4, gas entropy 
is significantly higher, and the cooling time is relatively 
long, even at very small radii. This suggests that while 
MKW 4 has a cool core, and may be developing towards 
a cooling flow, the core of AWM 4 has likely been heated. 
The gas fraction profile suggests that the inner halo has 
been 'puffed up' by the injection of energy, which is also 
responsible for the isothermal temperature profile. 

3. Discussion and Conclusions 

MKW 4 is similar to a number of systems observed by 
Chandra and XMM in that although its halo contains gas 
at a range of temperatures, there is no evidence o f cool- 
ing below ~0.5 keV. NGC 4636 lIXu et al.ll2002ft. M87 
ijSakelliou et alJ l2002j) and NGC 5044 l)Tamura et all 
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Fig. 5. — Three dimensional properties of the two clusters, as determined from their surface brightness and temperature profiles. AWM 4 
is shown in blue, MKW 4 in red. 



2003) all show this behaviour, with minimum temper- 
atures of 0.5-0.6 keV. Given this minimum temperature 
in MKW 4, it seems likely that the halo has been heated 
in the past and is now in the process of cooling and re- 
establishing a cooling flow. This heating could have been 



caused by some dynamical event, such as a sub-cluster 
merger, or by AGN activity. 

We calculated the time required for an isothermal halo 
to cool to the point where its temperature profile would 
match that which we observe in MKW 4. In the clus- 
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ter core, this cooling timescale is ^200 Myr. This is 
quite comparable to the activity cycle of an AGN, but 
very dissimilar to the spectroscopic age of the dominant 
galaxy. Along with the relaxed nature of the cluster, 
this rules out a subcluster merger as the source of heat- 
ing. We therefore suggest that MKW 4 is likely to have 
been heated by an AGN which is currently quiescent. We 
also note that the temperature profile of MKW 4 is not 
well described by either the 'universal' tem perature pro- 
file for relaxed clusters of lAllen et alJ l)200lf) or by models 
in wh ich cooling is balanced by conduction l)Voigt et al.l 
2002). This may support our suggestion that MKW 4 is 
not in a stable state and is likely to cool further. 

AGN heating is also likely responsible for the isother- 
mality of AWM 4. Taking the unlikely possibility that 
the isothermal temperature profile we observe is a sta- 
ble, long term feature of the system, an energy input of 
at least ~10 43 erg s _1 is required to balance the X-ray 
emission from the central 100 kpc. If, as is more likely, 
AWM 4 once had a temperature profile much like that 
of MKW 4, ~9xl0 58 erg would be needed to heat it to 
its current state. This is equivalent to an AGN inject- 
ing ^3xl0 43 erg s _1 for 100 Myr. We therefore believe 
that even considering efficiency factors, it is realistic to 
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assume AWM 4 has been heated by AGN activity. 

In conclusion, we find that AWM 4 and MKW 4 are 
very similar systems with two major differences: 

1. AGN activity: While AWM 4 is currently being 
heated by the AGN in NGC 6051, there is no on- 
going AGN activity in NGC 4073, and gas in the 
core of MKW 4 has been able to cool. 

2. Galaxy mass: NGC 6051 is considerably less 
massive than NGC 4073, leading to a less steeply 
peaked mass profile in the core of AWM 4. 

The observed differences in temperature and surface 
brightness profile between the two systems can be ex- 
plained largely through the effects of their respective cen- 
tral galaxies. The mass of these galaxies and the current 
phase of their AGN duty cycles appear to have profound 
effects on the structure and future evolution of the two 
clusters. 
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